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Abstract 



In this paper, we continue to study the blowup problem of the JV-dimensional compressible 
Euler or Euler-Poisson equations with repulsive forces, in radial symmetry. In details, we 
extend the recent result of "M.W. Yuen, Blowup for the Euler and Euler-Poisson Equations 
with Repulsive Forces, Nonlinear Analysis Series A: Theory, Methods & Applications 74 (2011), 
1465-1470." . We could further apply the integration method to obtain the more general results 
which the non-trivial classical solutions (p,V), with compact support in [0, R], where R > 
is a positive constant with p(t, r) = and V(t, r) = for r > R, under the initial condition 



where an arbitrary constant n > 0, blow up on or before the finite time T = 2R n+2 / (n(n + 
l)Ho) for pressureless fluids or 7 > 1. The results obtained here fully cover the previous known 
case for n = 1. 
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1 Introduction 



The compressible isentropic Euler (<5 = 0) or Euler-Poisson (6 — ±1) equations can be written in 
the following form: 

Pt+V ■ (pu) =0 

p[u t + (u ■ V)w]+VP =pV$ ( 2 ) 
A<f>(t,x) = Sa{N)p 

where a(N) is a constant related to the unit ball in R N . As usual, p = p(t, x) > and u = u(t, x) g 
H N are the density and the velocity respectively. P = P(p) is the pressure function. The 7-law 
for the pressure term P{p) could be applied: 



P(p) = Kpi 



(3) 



which the constant 7 > 1. If K > 0, we call the system with pressure; if K = 0, we call it 
pressureless. 

When 8 = — 1, the system is self-attractive. The equations ([2]) are the Newtonian descriptions of 
a galaxy in astrophysics [3] and [B]. When 6 = 1, the system is the compressible Euler-Poisson 
equations with repulsive forces. It can be used as a semiconductor model [8 . For the compressible 
Euler equation with 6 = 0, it is a standard model in fluid mechanics [PT] . And the Poisson equation 
([2]) 3 could be solved by 

(4) 



®(t,x) =<5 / G(x-y)p(t,y)dy 
Jr n 



where G is Green's function: 



G(x) 



\x\, N = l 
log |4 N ~2 
- 1 N > 3. 



(5) 



For the construction of analytical solutions for the systems, interested readers could refer to 



[15] . [18] . [PT] . [T6] and [24]. The local existence for the systems can be found in jT7]. [19], [2] and 
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P] . The analysis of stabilities for the systems may be referred to [23], [J, [T2J, [20], [2T], [22], 
[IH], El, [25], [SJ and @]. 

The solutions in radial symmetry could be in this form: 

p = p(t,r) and u = -V(t,r) =: -V (6) 
r r 

1/2 



with the radius r = f 

The Poisson equation ([2]) 3 becomes 



r" -1 ^ (t,a;) + (N- l)r N - z <S> r =a (iV) <5pr (7) 
^ = ^/%(M)^. (8) 



i) 



By standard computation, the systems in radial symmetry can be rewritten in the following form: 

N - 1 

Pt + Vp r + pV r H pV = 

(9) 

p(Vt + VV r ) + P r { P ) = P $r (/?)• 

In literature, Makino, Ukai and Kawashima first studied the blowup of tame solutions [20] for 
the compressible Euler equations (S = 0). Then Makino and Perthame investigated the corre- 
sponding tame solutions for the system with gravitational forces [21]. After that Perthame [22] 
obtained the blowup results for the 3-dimensional pressureless system with repulsive forces (S = 1). 
There are other blowup results for the systems in [12], [13], [5] and [4]. 

Very recently, Yuen |26] used the integration method to show that with the initial velocity 



H = [ rV dr>0, (Idi 
Jo 

the solutions with compact support to the Euler (5 = 0) or Euler-Poisson equations with repulsive 
forces (5 — 1) blow up in the finite time. 

In this article, we can observe that the condition (ITU|) could be more general to have the 
corresponding blowup results. In fact, we could further apply the integration method to extend 
Yuen's result as the following theorem: 

Theorem 1 Consider the Euler (5 = 0) or Euler-Poisson equations with repulsive forces (6 = 1) 
(0) in R . The non-trivial classical solutions (p, V) , in radial symmetry, with compact support in 
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[0, R], where R > is a positive constant (p{t,r) — and Vit, r) — for r > R) and the initial 
velocity: 

H = [ r n V dr > (11) 
Jo 

with an arbitrary constant n > 0, 

blow up on or before the finite time T — 2R n+2 /(n(n + l)Ho), for pressureless fluids (K = 0) or 
7 > 1. 

We remark that the condition 

p(t, r) = and V(t, r) = for r > R (12) 
in the theorem is called non-slip boundary condition [9] and [7]. 

2 Integration Method 

We just follow the integration method which was designed in [26] to obtain the further results. 
Proof. The density function p{t, x(t; x)) preserves its non-negative nature as we can integrate the 
mass equation <(2]) i : 

Dp , s 

-^+ P V-u = (13) 

with the material derivative, 

^=! +( - v > <i4 > 

to have: 

p{t, x) = p o (a;o(0, x Q )) exp f- V • u(t, a?(t; 0, x ))dt\ > (15) 

for poixoiO, Xq)) > 0, along the characteristic curve. 

Then we can manipulate the momentum equation @2 for the non-trivial solutions in radial 
symmetry, po ^ 0, to obtain: 

V t + VV r + K 1P ~<- 2 p r = $ r (16) 
V t + ^- r i\v 2 ) + K 1 pi-*p r = ® r (17) 
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r n V t + r n — (-V 2 ) + K-/r n p~<- 2 p r = r n <P r 
or 2 

with multiplying the more general function r™ with n > 0, on the both sides. 

We notice that this is the critical step in this paper to extend the blowup result in 

We could take the integration with respect to r, to equation (fT8|) . for 7 > 1 or K > 0: 

f R f R d 1 f R r R 

/ r n V t dr+ / r n — (-V 2 ) + / K ir n p^ 2 p r dr = / r n $ r dr 
Jo J dr 2 Jo Jo 



p(t, s)s"- l ds 



dr 



R r n V t dr+ f R r n ±{)-V 2 )+ ^ ^dp^ 1 > 
a r 2 Jo 7-1 

for (5 > 0. 

Then, the below equation can be showed by integration by part: 

J R r n V t dr - \ J Q R nr n - l V 2 dr + \ [R n V 2 {t, R) - 0" • V 2 {t, 0)] 
" Io R ^p^P^dr + [B«pi-\t, R) - 0" • pi-\t, 0)] > 0. 
The above inequality with the boundary condition (V(t, R) = and p(t, R) — 0), becomes 



R 1 f R f R K'ynr"' -1 

r n V t dr - - / nr n - x V 2 dr - / — '- p<~ x dr > 

2 Jo Jo I- 1 

- [ R r n Vdr~- I* n^-Wdr - f A ^" r "' 1 p i-^ dr > 
dtJo 2J io 7-1 

[ R Vdr^- 1 - l R -JL- V Hr^> f RK ^ 1 p^dr>0 
dtn+lj 2j a (n + l)r " J 7-I 

for n > and 7 > 1 or K = 0. 

For non-trivial initial density functions po > 0, we obtain: 

d 1 /■* 1 f R n 



Vdr n+1 - ± / ; V 2 dr n+1 > 



dtn + l J 2 Jo (n + l)r 

' /fl Fdr" +1 > ( R —^—V 2 dr n+1 > ( R V 2 dr n+1 



dtn+l Jo ~ Jo 2(n + l)r ~ 2(n + J 

- [ R Vdr n+1 > — [ R V 2 d 
dt Jo ~ 2R Jo 

We denote 



F := 



Jo n+lJo 



() 
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and apply the Cauchy-Schwarz inequality 

r-R 



V-ldr 



n+l 



r \ 1/2 { r R \ 1/2 

I / V 2 dr n+1 \ I J ldr n+1 \ 









ldr n+1 




J R V 2 dr n+1 ^j 









1/2 



n+l\ V 2 



f?Vdr<+ 1 



R- 



(R n+1 ) 

1/2 



I J V z dr 



rr2 



n(n + 1) H 2 n 
2i?"+ 2 ~2Rj Q 



for driving equation (|2"51) to be 



with inequality (f34]> . 



di V ' ~ 2RJ - 2R n + 2 



d H > n<y71 + 
di ~ 2R n + 2 



H > 



-2R n+2 H 



n(n + l)H t - 2R n + 2 ' 
Finally, we could require the initial condition 



H = / r n V o dr>0 
Jo 



(30) 
(31) 

(32) 

(33) 
(34) 



(35) 

(36) 
(37) 



(38) 



for showing that the solutions blow up on or before the finite time T — 2R n+2 / (n(n + I) Ho). 
This completes the proof. ■ 

We notice that the results in this paper fully cover the previous case for n = 1 [26 . Further 
researches are needed to have the corresponding results for the non-radial symmetric cases. 
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